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Summary 

This paper describes an integrated aerodynamic/ 
dynamic/structural (IADS) optimization procedure 
for helicopter rotor blades. The procedure combines 
performance, dynamics, and structural analyses with 
a general-purpose optimizer using multilevel decom- 
position techniques. At the upper level, the structure 
is defined in terms of global quantities (stiffnesses, 
masses, and average strains). At the lower level, the 
structure is defined in terms of local quantities (de- 
tailed dimensions of the blade structure and stresses). 

The upper level objective function is a linear 
combination of performance and dynamic measures. 
Upper level design variables include pretwist, point 
of taper initiation, taper ratio, root chord, blade 
stiffnesses, tuning masses, and tuning mass locations. 
Upper level constraints consist of limits on power 
required in hover, forward flight, and maneuver; 
airfoil drag; minimum tip chord; trim; blade natural 
frequencies; autorot at ional inertia; blade weight; and 
average strains. 

The lower level sizes the internal blade structure 
at several radial locations along the blade. The lower 
level optimization assures that a structure can be 
sized to provide the stiffnesses required by the up- 
per level and also assures the structural integrity of 
the blade. The lower level design variables are the 
box beam wall thicknesses and several lumped areas 
that are analogous to longitudinal stringers in a wing 
box cross section. The lower level objective func- 
tion is a measure of the difference between the upper 
level stiffnesses and the stiffnesses computed from the 
wall thicknesses and lumped areas. Lower level con- 
straints are on the Von Mises stress at the box corners 
for multiple-load cases generated by several flight 
conditions, limits on wall thicknesses for thin-wall 
theory, and other dimensional considerations. 

The IADS procedure provides an optimization 
technique that is compatible with industrial design 
practices in which the aerodynamic and dynamic 
designs are performed at a global level and the 
structural design is carried out at a detailed level 
with considerable dialogue and compromise among 
the aerodynamic, dynamic, and structural groups. 
The IADS procedure is demonstrated for several 
cases. 

Introduction 

Over the last decade, optimization techniques 
have been studied for application to the rotor blade 
design process. In reference 1, Miura presents a 
survey on the application of numerical optimization 


methods Jo helicopter design problems including ro- 
tor blade design. Most optimization procedures have 
dealt with a single discipline such as aerodynamics 
(refs. 2-4), structures (ref. 5), or dynamics (refs. 2 
and 6-9). However, the rotor blade design process 
is multidisciplinary involving couplings and inter- 
actions between several disciplines such as aero- 
dynamics, dynamics, structures, and acoustics. 
These couplings and interactions can be exploited 
by the optimization procedure if all the disciplines 
are accounted for simultaneously rather than sequen- 
tially. For instance, in a review (ref. 10) on the im- 
pact of structural optimization on vibration reduc- 
tion, Friedmann emphasizes the need to include the 
multidisciplinary couplings between aerodynamics, 
dynamics, and structures even when optimizing only 
for minimum vibration. 

Techniques and strategies for merging disciplines 
to obtain integrated rotorcraft optimization proce- 
dures are developing. In references 11 and 12, 
a plan is described for integrating the disciplines 
of aerodynamics, dynamics, structures, and acous- 
tics. As part of that plan, aerodynamics and 
dynamics have been incorporated systematically into 
performance (refs. 3 and 4) and airload/dynamic 
(ref. 13) optimization procedures resulting in an 
integrated aerodynamic/dynamic optimization pro- 
cedure (ref. 14). Reference 15 summarizes recent 
accomplishments based on that plan. 

Other multidisciplinary rotor blade optimization - 
work is described in references 16-19. References 16 
and 17 describe the formulation of a multidisciplinary 
approach to rotor blade design for improved perfor- 
mance and reduced fuselage vibrations. Reference 18 
describes a staged optimization procedure for a ro- 
tor for combined aerodynamics, dynamics, and struc- 
tures. Reference 19 describes a multidisciplinary 
optimization procedure to design high-speed prop 
rotors. 

What is lacking in previous multidisciplinary ro- 
tor blade optimization procedures is an efficient 
method to integrate structures or structural prop- 
erties. Usually, structures or structural properties 
are included in one of two ways — either as local de- 
sign variables (indirectly affecting the response of the 
blade) or as global design variables (directly affect- 
ing the response of the blade). When local design 
variables are used, the detail dimensions of a struc- 
tural member at one or more radial locations along 
the blade are used to generate structural properties. 
When global design variables are used, structural 
properties are the design variables. Both types of 
design variables have limitations. Using local design 
variables (e.g., refs. 6, 7, 18, and 19), such as wall 



thicknesses of the structural member, can lead to a 
large number of design variables that can be compu- 
tationally expensive. Also, this choice of design vari- 
ables is at odds with the traditional design practice in 
which chord, stiffness, and mass distributions along 
the blade are determined and then a structure is de- 
signed that matches these distributions. Using global 
design variables (e.g., refs. 2, 9, 13, 14, 16, and 17), 
such as stiffness and mass properties, in optimiza- 
tion also has disadvantages. When flapwise bending 
stiffness, chordwise bending stiffness, torsional stiff- 
ness, and extensional stiffness distributions are used 
as design variables, they are treated as independent 
quantities. In reality, these stiffnesses are not inde- 
pendent, and no guarantee can be given that a set 
of wall thicknesses can be found that will simultane- 
ously give these stiffnesses. 

This paper presents the methodology for incor- 
porating aerodynamics, dynamics, and structures in 
an integrated optimization procedure using both lo- 
cal and global design variables. Multilevel decom- 
position techniques based on reference 20 are used 
to add structural design variables and constraints 
to an existing aerodynamic/dynamic optimization 
procedure (ref. 14). The product is an integrated 
aerodynamic/dynamic/structural (IADS) optimiza- 
tion procedure. The multilevel formulation used in 
this paper was presented first in reference 15. An- 
other preliminary study of multilevel techniques ap- 
plied to rotor blade design is described in 
reference 21. 

The multilevel decomposition approach has been 
successfully applied to multidisciplinary problems 
(e.g., refs. 22-24). As originally proposed in refer- 
ence 25, the coordination procedure consisted of an 
optimum sensitivity analysis (ref. 26) and a set of 
equality constraints that relate the detailed (local) 
design variables of one subsystem to the global de- 
sign variables on the level above. However, as pointed 
out in reference 27, these equality constraints have 
caused difficulties in implementing multilevel decom- 
position procedures. The IADS procedure is based 
on the multilevel decomposition approach of refer- 
ence 20 which eliminates the equality constraints in 
the coordination procedure, thus allowing the use 
of the optimum sensitivity derivative found in refer- 
ence 28 that is less computationally costly. However 
in the IADS procedure, the set of lower level con- 
straints is replaced by an envelope function known as 
the Kreisselmeir-Steinhauser function (KS function, 
ref. 29) which further reduces the computational cost. 

First, the general multilevel decomposition strat- 
egy with two levels will be discussed. (Note that the 
systems with more levels are discussed in refs. 20, 22, 


and 25.) Next, the general strategy will be related to 
rotor blade design. Then, the IADS development in- 
cluding flowcharts of the upper and lower levels and 
the optimization procedure will be explained. Re- 
sults will be presented for several cases that demon- 
strate the strengths of the IADS procedure. 


Symbols and Abbreviations 

A area, ft 2 

AI autorotational inertia, 

£ Wjtj, lbm-ft 2 
j = l 


b 

C D 

Cl 

CF 

c d 

c d,all 

<Z,max 

q 

Cr 

Ct 

m p 

E 

EA 

El XX 

EI ZZ 

F 

fb,i 

fk 

fk,l 

fk,u 

ft,i 

A/ 

G 

GJ 


zth lumped area, ft 2 

box width, ft 

rotor coefficient of drag 

rotor coefficient of lift 

centrifugal force, lb 

airfoil section drag coefficient 

maximum allowable section drag 
coefficient 

largest section drag coefficient at 
azimuth angle 

airfoil section lift coefficient 

root chord, ft 

tip chord, ft 

pth upper level design variable 
Young’s modulus of elasticity, lb/ft 2 
extensional stiffness, lb 
chordwise bending stiffness, lb-ft 2 
flapwise bending stiffness, lb-ft 2 
lower level objective function 
zth bending frequency, per rev 
fcth frequency, per rev 

lower bound on A;th frequency, 
per rev 

upper bound on fcth frequency, 
per rev 

zth torsional frequency, per rev 

increment used in frequency win- 
dow, per rev 

torsional modulus of elasticity, lb/ft 2 
torsional stiffness, lb-ft 2 


2 


9c,i 

2 th lower level constraint function 

*^4,ref 

reference value of 4/rev rotating 

9i 

zth upper level constraint function 

vertical hub shear in forward 
flight, lbf 

9 max 

maximum lower level constraint 
function, maxj^} 

h 

fcth wall thickness, ft 

h 

t max 

nondimensional location of maxi- 

box height, ft 


mum airfoil thickness 

ITER 

number of trim iterations 

V(<j,t) 

Von Mises stress, lb/ft 2 

I XX 

chordwise moment of inertia, ft 4 

Vi 

2 th lower level design variable 

hz 

flapwise moment of inertia, ft 4 

W 

total blade weight, lbm 

J 

polar moment of inertia, ft 4 

Wj 

total weight of /th structural 
segment, lbm 

KS 

Kreisselmeir-Steinhauser function 

Xl 

nondimensional distance from airfoil 

ki 

2 th weighting factor in objective 
function 

3/7 r 

leading edge to left of wing box 
nondimensional distance from airfoil 

h 

factor of safety 


trailing edge to right of wing box 

MRA 

total number of aerodynamic 

Vi 

location of zth tuning mass 


segments 

Vtr 

point of taper initiation 

m t 

N 

2 th segment tuning mass, slug/ft 
number of blades 

Z l 

nondimensional lower airfoil 
coordinate 

NDV 

number of upper level design 
variables 

Z U 

nondimensional upper airfoil 
coordinate 


integer 

£ 

coordination parameter 

n 

£a 

number of constraint components in 

allowable average strain 

nc 



lower level 

£ y 

average strain 

nP 

frequency at n times the rotational 

$tw 

maximum pretwist, deg 


speed of the blade 

A 

Lagrange multiplier 

ns 

number of structural segments 

P 

pull-down factor 

OBJ 

upper level objective function 

a 

bending stress, lb/ft 2 

P 

main rotor power, hp 

°a 

allowable stress, lb/ft 2 

R 

blade radius from center of 

r 

shear stress, lb/ft 2 


rotation, ft 

V 

azimuth angle, zero over tail, deg 

T 

distance along blade from center of 

n 

rotor speed, rpm 


rotation, ft 

Subscripts: 

r j 

distance from center of rotation to 

a 

available 


center of j th segment, ft 

ff 

forward flight 

s N,ft 

N / rev rotating vertical hub shear in 
forward .flight, lbf 

h 

hover 



m 

maneuver 

S N, ref 

reference value of N / rev rotating 
vertical hub shear in forward 

max 

maximum 


flight, lbf 

min 

minimum 

^4,fiF 

4/rev rotating vertical hub shear in 

0 

optimum 

forward flight, lbf 

ref 

reference 
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Superscripts: 

i ith component 

L lower level 

T transpose 

U upper level 

A bar over a symbol indicates a nondimensional 
quantity, and an asterisk (*) used as a superscript 
indicates an upper level design variable. 

Multilevel Optimization Strategy 

With a multilevel decomposition approach 
(refs. 20, 22, and 25), a large complex optimiza- 
tion problem is broken into a hierarchy of smaller 
optimization subproblems. This hierarchy can be 
thought of as levels of increasing detail. At the upper 
level, the subproblem is formulated in terms of global 
quantities that describe the overall behavior of the 
entire system. On the lower level, the subproblems 
are stated in terms of local quantities and local con- 
straints that have only a small impact on the entire 
system. Each of these subproblems use local design 
variables to reduce the violation of constraints that 
are unique to that subproblem. The coupling be- 
tween the upper level subproblem and the lower level 
subproblems is preserved through a coordination pro- 
cedure such as that described in references 20 or 25. 
This coupling represents a dialogue between the lev- 
els that, upon convergence, establishes compatibility 
between the two levels. 

Figure 1 illustrates a generic two-level optimiza- 
tion procedure. Note that the analysis proceeds from 
the upper level to the lower level while the optimiza- 
tion proceeds from the lower level to the upper level. 
First, the upper level analysis initializes all the global 
quantities and responses and then provides informa- 
tion to each lower level subproblem. Next, individual 
lower level optimizations are performed that reduce 
local constraints as much as possible and that provide 
information to the coordination procedure. Finally, 
the upper level optimization occurs. The preceding 
description defines 1 cycle. This entire process is re- 
peated for several cycles. Convergence occurs when 
all the constraints (both upper level and lower level) 
are satisfied and the upper level objective function is 
minimized. 

The rotor blade optimization problem can be de- 
composed into one subproblem affecting the global 
response of the blade and three subproblems affect- 
ing portions of the blade. Quantities such as power 
required, blade trim, autorotational inertia, natural 


frequencies, total blade weight, and average strain 
describe the global response of the blade. The entire 
blade must be analyzed to obtain these response 
quantities. Quantities such as stresses are detailed 
response quantities since only a portion of the blade 
must be considered to obtain these response quanti- 
ties. Therefore, a two-level decomposed rotor blade 
optimization problem can be defined as shown in fig- 
ure 2. The upper level optimizes the blade by chang- 
ing global quantities such as blade planform, twist, 
and distributions of mass and stiffness. The upper 
level chord, mass, and stiffness distributions are 
treated as independent quantities. The reconciliation 
between these distributions is done on the lower level, 
which consists of several independent subproblems at 
stations along the blade radius. These subproblems 
optimize detailed cross-sectional dimensions to sat- 
isfy stress constraints and to reconcile the upper level 
independent mass, chord, and stiffness distributions 
with the lower level calculated mass and stiffness dis- 
tributions. This reconciliation is improved further by 
a set of upper level coordination constraints. (See 
appendix A.) First, the upper level analysis and op- 
timization will be described, then the lower level 
analysis and optimization, and last the overall IADS 
system. 

Upper Level Analysis and Optimization 

The purpose of the upper level analysis is to 
evaluate the overall rotor blade design on the ba- 
sis of performance, dynamic, and global structural 
measures. (For a description of the rotor blade de- 
sign philosophy, see refs. 3, 4, 11, 12, 14, and 15.) 
The upper level analysis is similar to the integrated 
aerodynamic/dynamic analysis reported in refer- 
ence 14 with the addition of extensional stiffness de- 
sign variables, strain constraints, and coordination 
constraints. As shown in figure 3, the blade is evalu- 
ated for three flight conditions: hover, forward flight, 
and maneuver. The Langley-developed hover anal- 
ysis program HOVT (a blade-element momentum 
analysis based on ref. 30) is used to predict power re- 
quired in hover. The comprehensive helicopter anal- 
ysis program CAMRAD/JA (ref. 31) is used to pre- 
dict rotor performance (e.g., trim, airfoil drag, and 
power required), loads, and frequencies for forward 
flight and maneuver. The maneuver flight condition 
simulates a coordinated turn in terms of an increased 
load on the forward-flight lift requirement. 

The rotor blade design process is defined in terms 
of aerodynamic performance, dynamics, and global 
structural requirements. Satisfactory aerodynamic 
performance is defined by the following four require- 
ments. First, the power required for any flight 
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condition must be less than the available power. 
Second, airfoil section drag along the blade radius on 
the advancing and retreating side of the rotor disk in 
both forward flight and maneuver must be less than 
a maximum allowable value. Third, the rotor must 
trim at each flight condition. The rotor is trimmed 
to a constant lift in forward flight and a (different) 
constant lift in maneuver which ensures that the ro- 
tor has no loss in lift capability or maneuverability 
even if solidity decreases from the initial to the fi- 
nal design. Incorporation of a maneuver flight con- 
dition is used in place of a constraint on solidity, be- 
cause low-speed maneuver determines rotor solidity 
(ref. 32). Fourth, the blade tip chord must be larger 
than a prescribed minimum value. Satisfactory dy- 
namics is defined in terms of limits on vibrational 
frequencies. The blade is designed so that the natu- 
ral frequencies (both bending and torsional) do not 
coincide with integer multiples of the rotor speed. 
Also, the blade must have sufficient autorotational 
inertia as a safety measure needed in case of engine 
failure. In addition to satisfying these design require- 
ments, the blade weight must not exceed some upper 
limit. Satisfactory structural requirements are de- 
fined in terms of limits on the average axial strains for 
forward flight and maneuver flight conditions. The 
upper level optimization problem is formulated next 
in terms of design variables, objective function, and 
constraints. 

Upper Level Design Variables 

The upper level design variables are the blade 
planform, stiffnesses, and tuning masses. (See fig. 4.) 
The blade planform is defined by the point of taper 
initiation ( ytr ), root chord (cr), taper ratio (cv/q), 
and maximum pretwist (0 tw ). The blade is rectan- 
gular from the root to y tT and then tapers linearly 
to the tip. The pretwist varies linearly from the 
center of rotation to the tip. Global design vari- 
ables include the blade chord wise, flapwise, torsional, 
and extensional stiffnesses (denoted by EI XX , EI ZZ , 
GJ, and EA , respectively) at three radial locations: 
blade root, point of taper initiation, and blade tip. 
The stiffnesses are assumed to vary linearly between 
these points and are treated as independent quanti- 
ties. The remaining design variables are three tun- 
ing masses (denoted by mi, m2, and m3) and their 
locations (denoted by y\,y2, and 2/3), respectively. 
The total blade mass consists of the structural mass 
(which is assumed to be constant) plus the sum of 
the tuning masses. No attempt is made to reconcile 
the change in weight with the change in design vari- 
ables because the present work is based on extending 
the procedure of reference 14 to include structures. 
However, this reconciliation is possible. (See ref. 15.) 


The center of gravity and aerodynamic offsets are co- 
incident with the blade elastic axis. The number of 
blades, rotor radius, rotational velocity, airfoils, and 
airfoil distribution are preselected and fixed. 


Upper Level Objective Function 

The objective function to be minimized is a com- 
bination of performance and dynamics measures and 
is formulated as 


OBJ = kll p^+k 2 ^+k 3 -^+k 4 Sn * 


Ph, ref Z Pff,ref 


m, ref 


S N t re f 


( 1 ) 


where P/>,P ff, and P m are the powers required in 
hover, forward flight, and maneuver, respectively. 
The symbol N is the number of blades, and 
is the N / rev rotating vertical hub shear in forward 
flight. The terms &i, ^2? and are weighting fac- 
tors chosen by the user, and P hM , P ff ref , P m ref , 
and S jv,ref are reference values used to normalize 
and nondimensionalize the objective function com- 
ponents. The usefulness of this objective function 
was demonstrated in reference 14. 


Upper Level Constraints 

The upper level constraints are grouped into 
performance, dynamic, structural, and coordination 
constraints. This section of the paper discusses the 
performance, dynamic, and structural constraints. 
The coordination constraints are discussed later in 
the paper. The performance and dynamic constraints 
are the same as those used in reference 14. By con- 
vention, the ith constraint gi is satisfied if it is less 
than or equal to zero. 

Performance constraints . The performance 
constraints are on power required, trim, airfoil sec- 
tion drag, and blade tip chord. The requirement that 
the power required be less than the power available 
is given by 

9i = ^ - 1 < 0 (2) 

r a 

for each flight condition, where Pj is the power 
required for the zth flight condition and P a is the 
power available. 

The requirement on the airfoil section drag trans- 
lates into a constraint that each airfoil section distrib- 
uted along the rotor blade operate at a section drag 
coefficient less than a specified allowable value 
c d,all- (See appendix B.) This leads to 24 constraints 
per flight condition because the blade is analyzed in 
azimuth increments of 15° around the rotor disk. At 
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(6a) 


a given azimuth angle ('!'), the constraint is formu- 
lated as 

— l < 0 ($ = 15 °, 30 °, 45 °,... , 360 °) (3) 

Cd, all 


an amount A/. Thus, for the upper bound, 

9i = -p- - 1 < 0 

Jk,u 

and for the lower bound, 


where c^aii is the allowable drag coefficient and 
c dmax * s l ar S es f drag coefficient at any radial 
station. (Note that the drag coefficients in the 
reverse-flow region occurring on the retreating side 
of the rotor disk are ignored.) In the present work, 
the same value for c^&w is used on the advancing and 
retreating side of the rotor disk. This simplifying 
assumption can easily be lifted. 

The trim requirement is difficult to translate into 
a mathematical constraint. The trim constraints 
in forward flight and maneuver are implemented by 
using the method developed in reference 3, which 
expresses the constraint in terms of the number 
of trim iterations (ITER), the maximum number 
of trim iterations allowed (ITE R ma x)^ and the pth 
nondimensional design variable (DV P ). The heuristic 
trim constraint is given by 


9i - (ITER - ITER max + 1) 



<0 (4) 


where NDV is the number of design variables. In 
the development of this equation in reference 3, the 
addition of the summation term was found to im- 
prove convergence because it allowed calculation of 
the change in the trim constraint with respect to the 
change in a single design variable. 

The final performance requirement is a constraint 
used to ensure that the blade tip chord does not 
become too small. Thus, 

gi — 1 — £ 0 (5) 


9i = l ~ 


fk 

~~ < 0 
JkJ 


(6b) 


where /*. u has a value that is A / below n+ 1 per rev 
and has a value that is A / above nj rev for the 
applicable n. For example, if A / is 0.1/rev and is 
5.6/rev, then nP would be 5/rev and (n + 1 )P would 
be 6/rev. Thus, f 4 jU and would be 5.9/rev and 
5.1/rev, respectively. Formulating the constraints in 
this manner allows the frequencies to change from 
one optimization cycle to the next cycle provided 
the frequencies avoid approaching integer multiples 
of the rotor speed. This formulation is different from 
the approaches used in references 13, 16, and 17 
in which the frequencies are kept within prescribed 
windows based on the reference blade frequencies. 
In this work, constraints are placed on frequencies 
in both forward flight and maneuver because blade 
collective pitch and the amount of modal coupling 
may be different for the two flight conditions, and 
therefore the frequencies can be different. 

The constraint that the blade weight should be 
less than some maximum value is formulated as 


9i = 


W 

Umax 


- 1 < 0 


( 7 ) 


where W is the total blade weight and W mixx is the 
maximum allowable weight. The total blade weight is 
the structural mass distribution (which is constant) 
plus the sum of the tuning masses. 

Finally, the blade must have enough auto- 
rotational inertia (AI) for safe autorotation in case 
of engine failure. The constraint is formulated so 
that the autorotational inertia of the blade is greater 
than some minimum value AI m } n . Thus, 


where Ct is the tip chord and c $ |in m is the minimum 
tip chord allowed. This is a practical constraint used 
to assure validity of the airfoil tables and to address 
manufacturing considerations. 

Dynamic constraints . The dynamic con- 
straints are on frequencies, total blade weight, and 
autorotational inertia. The constraint on the fcth 
frequency (either a bending or a torsional fre- 
quency) is formulated such that the frequency is sep- 
arated from integer multiples of the rotor speed by 


9i = 1 - 


AI 

Almii 


< 0 


( 8 ) 


Structural constraints . The structural con- 
straints are on the average axial strains. The struc- 
tural constraints evaluated at the same radial loca- 
tions that are used to define the design variables 
(fig. 4) are imposed on the average axial strains (e y ) 
as follows: _ 

Pi = — — 1 < 0 (9a) 

£a 


6 


and 


9i = - 1 - ^ < 0 (9b) 

£a 

where e a is the magnitude of the allowable strain and 


l f C F 
Sy ~ EA 


( 10 ) 


where CF is the centrifugal force, EA is the ex- 
tensional stiffness, and If is a safety factor on the 
loads. The strain constraints are calculated using 
loads from both the forward flight and the maneuver 
flight conditions. 


Upper Level Optimization 

The upper level optimization consists of the 
general-purpose optimization program CONMIN 
(ref. 33) and an approximate analysis used to re- 
duce the number of HOVT and CAMRAD/JA anal- 
yses during the iteration process. The approximate 
analysis is used to extrapolate the upper level objec- 
tive function and upper level constraints with linear 
Taylor-series expansions using derivatives of the ob- 
jective function and constraints with respect to the 
design variables 


OBJ = OBJ 0 + 


NDV T 

E oOBJ 

dDVi 

2=1 1 


ADV, 

o 


(ii) 


and 


9 = 90 + 


NDV 

E 


2=1 


dg 

dW 1 


ADVj 


( 12 ) 


The assumption of linearity is valid over a suitably 
small change in the design-variable values and will 
not introduce a large error into the analysis provided 
that the changes ADV are small. Errors that may 
be introduced by use of the approximate analysis 
are controlled by imposing “move limits” on each 
design variable during the iteration process. A move 
limit that is specified as a fractional change of each 
design-variable value is imposed as an upper and 
lower design- variable bound. At the present time, 
the move limits are manually adjusted. 


Lower Level Analysis and Optimization 

This section of the paper describes the lower level 
analysis and lower level optimization procedure. The 
purpose of each lower level optimization is to assess 
whether a structure at the given radial location can 
be sized to provide the stiffnesses required by the 
upper level optimization and still have the strength 
to withstand loads calculated by the upper level 


analysis. The lower level optimizations can be done 
in parallel because they are independent. 

For simplicity, because closed-form equations can 
be derived (see appendix C), the structural member 
(fig. 5) is assumed to be a thin-walled isotropic box. 
The box cross section is symmetric about the hori- 
zontal axis with wall thicknesses (£j) and lumped ar- 
eas (aj) which are analogous to longitudinal stringers 
in a wing box cross section. The outer dimensions b 
(the box beam width) and h (the box beam height) 
are functions of the upper level design variables be- 
cause b and h depend on the local chord and the 
local airfoil thickness. The values of b and h are de- 
termined by placing a box of maximum area within 
the airfoil cross section by using the method of ref- 
erence 34. (See appendix D.) 

Lower Level Design Variables 

The design variables are the three wall thicknesses 
and £ 4 ) and the three lumped areas (a\,a 2 , 
and < 23 ). The lumped areas are used to give the 
lower level more flexibility in matching the upper 
level stiffnesses. For the present implementation, the 
lumped areas are assumed to be square areas. 

Lower Level Objective Function 

The objective function is a measure of the dif- 
ference between the stiffnesses required on the up- 
per level and those determined from the lower design 
variables 


F = 


\EI Z . 


(EI ZZ ) 


* 1 2 


+ 


(Ei zz y 

gj - ( G jy 


+ 


El XX ( El XX ) 


: 1 2 


(Ei xx y 


0 gj y 


(13) 


where a starred quantity ( )* denotes an upper level 
design variable. The lower level cross-sectional prop- 
erties 1 xx, I zz , and J are computed (see appendix C), 
E is Young’s modulus of elasticity, and G is the tor- 
sional modulus of elasticity. 


Lower Level Constraints 

The constraints are enforced on the extensional 
stiffness, stresses, and the physical dimensions of the 
wall thicknesses and lumped areas. The extensional 
stiffness constraint that requires the lower level cal- 
culated extensional stiffness EA (appendix C) to be 
equal or greater than the upper level extensional stiff- 
ness (EA)* (an upper level design variable) is given 

by 

s » = 1 - m? * 0 < 14 > 
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at the given cross section. The extensional stiffness 
appears in a constraint rather than in the objective 
function (eq. (13)) where the other stiffnesses appear. 
This is done because the role of EA in the upper level 
is limited to satisfying the strain constraints (eq. (9)). 
The lower level is responsible only for assuring that 
the value of EA is at least as large as the value 
needed in the upper level; i.e. , close matching of EA 
to (EA)* is not required. 

The stress constraint evaluated at the corner of 
the box cross section shown in fig. 5 has the form 


9c,i ” 


Vfar) 


- 1 <0 


(15) 


where a is the bending stress, r is the shear stress, 
and V(a, r) is the Von Mises stress measure. (See 
appendix C.) Two stress constraints are used: in one, 
r is based on the vertical wall thickness, and in the 
other, r is based on the horizontal wall thickness. 

A set of constraints is imposed on the lower level 
wall thicknesses to assure that the section remains 
a thin-walled section and that the expression for J 
remains valid. (See appendix C.) These constraints 
are 


to 

-4r-l <0 
0.16 

(j = 2 and 4) 

(16) 

t-i 

\ 1 < 0 
O.lh 

(j = 1 and 3) 

(17) 


where b and h are the width and height of the box 
cross section, respectively. 

A set of constraints is imposed on the lumped ar- 
eas and wall thicknesses that require that the dimen- 
sions are physically possible (i.e., that the lumped 
areas can fit inside the box cross section). These 
constraints are 

9c, i = - ( *~ ~ V^T- §>/g) ^ 0 (18) 

(19) 

g c j = — (h — t\ - <3 — 2 v /o7) < 0 (20) 

9c,i = - {h - h ~ h - 2 ^/ 02 ) < 0 (21) 

g c j = — {h — t\ — £3 — 2 ^/ 03 ) < 0 (22) 

In addition, a set of constraints representing up- 

per and lower bounds on the design variables is used. 


For the kth design variable, the lower bound is given 
by 

9c,i = v k ,i -v k <0 (23) 

and the upper bound is given by 

9c, i =n~ v k,u < 0 ( 24 ) 

where v k \ and v kiU are the lower and upper design 
variable bounds, respectively. 

For convenience, the set of lower level constraints 
defined by equations (14)— (24) is replaced by a single 
cumulative constraint, an envelope function known 
as the KS function (ref. 29), which approximates the 
active constraint boundary 


KS — <?max H" l n 
P 


E e 

-i=l 


.p(9c t i~9m ax) 


< 0 (25) 


where g max is the maximum constraint component 
from equations (14)-(24), nc is the number of lower 
level constraint components, and p is defined by 
the user. Initially, p is small and then increases 
until a maximum value p m ax is reached. For large 
values of p, the value of KS approaches p m ax- The 
KS function is a single measure of the degree of 
constraint satisfaction or violation and is positive 
(violated) if at least one of the constraints g Ci i is 
violated. The KS function is a single-valued function 
that is continuous and differentiable. This property 
becomes important when implementing the upper 
and lower levels as described in the section on the 
overall organization of the IADS procedure. 


Lower Level Optimization Procedure 

The flowchart for each lower level optimization 
procedure is shown in figure 6. Loads, local chord, 
box beam width, box beam height, and upper level 
stiffnesses are passed down from the upper level anal- 
ysis. The lower level design variables (fig. 5) are 
used to calculate lower level stiffnesses. Von Mises 
stresses are calculated using the loads from the for- 
ward flight and maneuver analyses. The lower level 
objective function (eq. (13)) and cumulative con- 
straint (eq. (25)) are evaluated. The lower level op- 
timizations are performed using the general-purpose 
optimization program CONMIN. Exact analyses are 
used to evaluate the objective function, the con- 
straint, and any gradients computed by CONMIN. 
The optimization process is converged when the ob- 
jective function is minimized and the cumulative con- 
straint is satisfied. After convergence, the process 
returns to the upper level. 
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Coordination Between Upper and Lower 
Levels 

The coordination between upper and lower lev- 
els is implemented by upper level constraints. These 
constraints are imposed to encourage changes in the 
upper level design variables that promote consistency 
between the upper and lower level stiffnesses. Specif- 
ically, these constraints (one for each lower level op- 
timization) have the form 

g = F U -(l+e)F^ <0 (26) 

where is the most recent value of the lower 
level objective function (i.e., the optimum value of 
eq. (13)), F u is an estimate of the change in Fq that 
would be caused by a change in the upper level design 
variable values, and e is a specified tolerance defined 
as the coordination parameter. (See appendix A.) 
The importance of this parameter will be discussed 
later. 

Equation (26) is the general form of the coordi- 
nation constraint as formulated in reference 15. As 
shown in appendix A, the coordination constraint 
can be approximated in terms of the lower level to- 
tal optimum sensitivity derivative that expresses how 
the optimum lower level objective function and lower 
level active constraint will change with a change in 
upper level design variable. 

Overall Organization of IADS Procedure 

The conceptual IADS procedure is shown in fig- 
ure 2. It consists of an upper level analysis (fig. 3), 
three lower level optimizations (fig. 6), and a coor- 
dination task. The actual IADS procedure is more 
complicated and requires, in addition, an upper level 
sensitivity analysis and three lower level optimum 
sensitivity analyses. 

The flowchart for the IADS procedure is shown in 
figure 7. First, the upper level analysis is executed 
for the current set of design variables providing all 
the information needed to calculate the upper level 
objective function and constraints with the exception 
of the coordination constraints. The upper level 
analysis also provides the loads, local chord, box 
beam width, box beam height, and stiffnesses (to 
be matched) to the lower level analysis. Each lower 
level optimization is performed to obtain a set of 
lower level design variables that match the current 
upper level bending and torsional stiffnesses as close 
as possible. 

Next, an upper level sensitivity analysis is per- 
formed consisting of forward finite-difference deriva- 
tives (or gradients) of the upper level analysis. These 


derivatives are required to approximate the upper 
level objective function and upper level constraints 
during the upper level optimization. In addition, the 
loads and local chords corresponding to the changes 
in the upper design variables are saved. These quan- 
tities are used in the three lower level optimum sensi- 
tivity analyses to approximate the coordination con- 
straint (eq. (26)). Appendix A describes how the 
coordination constraint is expressed in terms of the 
total optimum sensitivity derivative involving both 
changes in the optimum lower level objective func- 
tion with respect to changes in the upper level de- 
sign variables and changes in the active lower level 
constraint with respect to changes in the upper level 
design variables. 

Finally, the upper level optimization occurs con- 
sisting of CONMIN and an approximate analysis. 
This describes 1 cycle of the IADS procedure. The 
process is repeated for additional cycles until conver- 
gence is achieved. A very strict convergence crite- 
rion is used for demonstration purposes. The overall 
procedure is converged when the change in the up- 
per level objective function is less than 0.5 x 10“ 5 
over three consecutive cycles and all the constraints 
(both upper and lower level) are satisfied. A step 
size of 0.001 is used to compute the finite-difference 
derivatives. 

Demonstration of IADS Procedure 

This section of the paper describes the analytical 
blade model, the mission definition, the optimization 
problem, and the optimization results used to demon- 
strate the IADS procedure. Results are presented for 
three studies: (1) the effect of initial design, (2) the 
effect of the coordination parameter e , and (3) the 
comparison between a single-level and multilevel op- 
timization approach. 

Analytical Blade Model 

The analytical blade model used to demonstrate 
the IADS procedure represents a wind tunnel model 
of a rotor blade for a four-bladed helicopter having 
a blade radius of 4.68 ft. Three sets of advanced ro- 
torcraft (RC) airfoils are used along the blade: the 
RC(4)-10 airfoil (ref. 35) from the root to 85 per- 
cent radius, the RC(3)-10 airfoil (ref. 36) from 85 to 
95 percent radius, and the RC(3)-08 airfoil (ref. 36) 
from 95 percent radius to the tip. Tables of exper- 
imental two-dimensional airfoil data for these three 
airfoil types are used in both HOVT and 
CAMRAD/JA. The analytical model of the blade 
uses 19 aerodynamic segments for HOVT, and it 
uses 50 structural segments and 18 aerodynamic seg- 
ments for CAMRAD/JA. HOVT is used to predict 
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the power required in hover using nonuniform inflow 
(no wake is included) by trimming to a constant lift 
(Cl). CAMRAD/JA is used to predict rotor perfor- 
mance, loads, and frequencies using uniform inflow 
with empirical inflow correction factors for the for- 
ward and maneuver flight conditions. Uniform in- 
flow is used to save on computational costs. (Note 
that even though an approximate analysis is used 
in the upper level optimization, 46 CAMRAD/JA 
analyses are required per optimization cycle.) In 
CAMRAD/JA an isolated rotor analysis is used that 
trims the rotor to constant lift (Cl) and drag (Cp) 
and zero flapping angle relative to the shaft using col- 
lective, lateral cyclic, and longitudinal cyclic pitch. 
From the modal analyses in CAMRAD/JA using 
10 bending modes and 5 torsional modes, only the 
first 6 bending frequencies are below 10 per rev and 
need to be constrained for a four-bladed rotor. Be- 
cause corresponds to a rigid-body mode and 
is the 1/rev frequency, the first two frequencies are 
not constrained. Constraints are placed on the first 
four bending frequencies (/& } 3, 4, and fog are flap- 
ping dominated and is lead-lag dominated) and 
the first two torsional frequencies (/^i represents the 
rigid-body torsional mode due to the control system 
stiffness and 2 represents the first elastic torsional 
mode). 

Mission Definition 

The flight conditions are a constant lift of 1# 
(331 lb and Cl = 0.0081), a propulsive force of 32 lb 
(C D = -0.000811), and an advance ratio of 0.35 for 
the forward flight condition and a constant lift of 
401 lb (Cl = 0.00985), a propulsive force of 23 lb 
(Cd - —0.000596), and an advance ratio of 0.3 for 
the maneuver flight condition, which is for a load 
factor of 1.22. These flight conditions and the load 
factor are similar to those used in reference 37. 

Optimization Problem 

The objective function is a combination of the 
power required in hover, forward flight, and maneu- 
ver and of the 4/rev rotating vertical hub shear in 
forward flight. The objective function is chosen to 
be one dominated by performance with little empha- 
sis on dynamics. Of the three powers, reducing the 
power required in hover is assumed to be the most im- 
portant; it will have twice the weight as the other two 
powers. Several values were tried for the weighting 
factor on the hub shear term. To obtain the proper 
balance between performance and dynamics, £4 must 
be between one and two orders of magnitude less 
than k v Thus, for this case, the weighting factors 


are chosen to be k\ — 10, k 2 = = 5, and k 4 = 0.5. 

Therefore, we have 

OBJ = 10-5^ + 5-^- + 5-^- + 0.5|^ (27) 

ref Mf,ref ^m^ef *^4, ref 

where /\ re f, -Pff,ref> ^m.reH anci 54, ref are 15 hp, 
13 hp, 12 hp, and 2 lbf, respectively. The reference 
values are chosen to be representative of the powers 
required and the hub shear for all the initial blade 
designs used in this work. 

The upper and lower bounds for the design vari- 
ables are given in table 1. On the upper level, 22 de- 
sign variables and 95 constraints are used. On the 
lower level, 6 design variables and 1 cumulative con- 
straint (the KS function with 24 components) are 
used at each of the 3 spanwise locations (i.e., the 
root, the point of taper initiation, and the tip). 

Parameters and flight conditions are summarized 
in table 2. Because the blade is made of aluminum, E 
has a value of 15.26 x 10 8 lb/ft 2 , the allowable strain 
s a has a value of 0.05 ft /ft, and the allowable stress 
a a is 8.352 x 10 6 lb/ft 2 . The values for minimum 
tip chord (c t]in j n ), power available ( P a ), minimum 
autorotational inertia, and maximum allowable drag 
coefficient (c^an) are 0.083 ft, 20 hp, 23.69 lbm-ft 2 , 
and 0.12, respectively. Frequencies must be at least 
0.1 away from a per-rev value (A/ = 0.1/rev in 
eq. (6)). 

Study on Effect of Initial Designs 

The IADS multilevel optimization procedure is 
demonstrated for three examples using the three 
starting points shown in figure 8. Example 1 
(fig. (8a)) uses a rectangular planform with a pretwist 
of -9°, a root chord of 0.3449 ft, and upper level stiff- 
ness design variables initialized to be consistent with 
the lower level initial wall thickness and lumped areas 
(i.e., matched stiffnesses). Example 2 (fig. (8b)) uses 
a tapered planform with a pretwist of —16°, a root 
chord of 0.45 ft, and matched stiffnesses. The blade 
is rectangular to 80 percent radius and then tapers 
linearly to the tip with a 3-to-l taper ratio. Exam- 
ple 3 (fig. (8c)) uses the same planform and pretwist 
as example 2 except that the upper and lower level 
stiffnesses are unmatched. All these examples use a 
value of -0.4 for the coordination parameter (e) in 
equation (26). The importance of the choice of e is 
examined in a later section of the paper. 

Example 1: rectangular planform (“ini- 

tially matched stiffnesses”). The starting point 
for the optimization is the rectangular blade shown 
in figure 8(a). The upper and lower level stiffnesses 
are matched because the upper level stiffnesses are 
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started with the stiffnesses determined by the ini- 
tial lower level design variables. This is an infeasi- 
ble starting point because the lower level stress con- 
straints at the root are violated. (See the results 
given in table 3.) The initial and final values for 
the blade planform, performance measures, and dy- 
namics measures are given in table 3(a). The initial 
and final values for the constrained frequencies are 
given in table 3(b). Notice that the final value for 
the fourth bending frequency f h A is in a different fre- 
quency range than the initial value. Final values for 
the lower level design variables and the upper level 
stiffnesses are given in table 3(c). The final design is 
able to improve the performance characteristics from 
the initial blade and satisfy all the constraints. Com- 
pared with the initial values, the final design repre- 
sents a 2.1-, 2.3-, 2.3-, 47.6-, and 3.2-percent reduc- 
tion in the power required in hover, forward flight, 
and maneuver; hub shear; and upper level objective 
function, respectively. 

The final stiffness distributions for the upper 
(required values) and lower levels (actual values) 
are shown in figure 9. The matching of the chord- 
wise bending stiffness (EI XX ) (fig. 9(a)), the flapwise 
bending stiffness ( EI ZZ ) (fig. 9(b)), and the torsional 
stiffness ( GJ ) (fig. 9(c)) are extremely good. As 
shown in figure 9(d), the lower level is able to ob- 
tain an extensional stiffness distribution higher than 
the minimum requirement set by the upper level. 

Convergence histories of the individual terms 
of the lower level objective function (eq. (13)) are 
shown in figure 10 for the three locations: the root 
(fig. 10(a)), the point of taper initiation (fig. 10(b)), 
and the tip (fig. 10(c)). Each term (denoted stiff- 
ness deviation) is a measure of how well the upper 
and lower stiffnesses match. Initially, the stiffnesses 
are matched, but the stress constraints are violated 
at the root. Therefore, the lower level design vari- 
ables must change to satisfy these constraints while 
keeping the upper and lower level stiffnesses matched 
as close as possible. Notice that_the chordwise stiff- 
ness at the root, the torsional stiffness at the point of 
taper initiation, and the flapwise stiffness at the tip 
are the last stiffnesses to match. Further, it appears 
that stiffnesses at the point of taper initiation are 
particularly difficult to match. This difficulty may 
be due to the fact that the point of taper initiation 
is a design variable but the root and tip positions are 
fixed. 

The reason for the deviations in the stiffness is 
that the upper and lower levels are in conflict. One 
component of the upper level objective function is 
the hub shear which can be reduced significantly 
by increasing the blade stiffnesses. On the upper 


level, if the optimizer did not have to be concerned 
with stiffness matching, it would increase the upper 
level stiffnesses. Without the lower level to keep the 
stiffnesses in check, a heavy or nonbuildable blade 
might result. 

The information shown in figure 10 is collected 
and used to determine when a move-limit adjustment 
is necessary for an upper level design variable during 
the overall optimization process. (Recall that an ap- 
proximate analysis is used on the upper level and an 
exact analysis is used on the lower level.) At the 
present time, no automatic move- limit adjustment 
is used in the approximate analysis on the upper 
level. Instead, the IADS procedure is run for 8 cy- 
cles and then the stiffness deviations are examined. 
When the stiffness deviation increases (e.g., cycle 16), 
the design- variable move limits are manually reduced 
and the optimization process is continued for another 
8 cycles. In practical applications, the optimiza- 
tion procedure would terminate after about 30 cy- 
cles; however, for demonstration purposes, the con- 
vergence criterion is set to a very small value. Both 
the upper and lower levels have the same tight con- 
vergence criterion on each cycle. Overall convergence 
of the IADS procedure might improve if the conver- 
gence criterion is relaxed initially and then tightened 
as the optimization proceeds. 

Example 2: tapered planform (“initially 

matched stiffnesses”). The starting point for 
the optimization is the tapered blade shown in fig- 
ure 8(b). Initially, the upper and lower level stiff- 
nesses are matched because the upper level stiffness is 
determined by the lower level design variables. How- 
ever, this is an infeasible starting point because a 
thin-wall-theory constraint is violated on the lower 
level. The initial and final values for the blade plan- 
form, performance measures, and dynamics measures 
are given in table 4(a). The initial and final con- 
strained frequencies are included in table 4(b). The 
final design is able to improve the performance char- 
acteristics from the initial blade. However, the hub 
shear increases from the initial value. 

Figure 11 shows the final stiffness distributions 
for the upper levels (required values) and lower levels 
(actual values) for the chordwise bending stiffness 
(fig. 11(a)), flapwise bending stiffness (fig. 11(b)), 
and torsional stiffness (fig. 11(c)). As shown in 
the figure, the stiffness matching is good, although 
not as good as in example 1. The lower level is 
able to obtain an extensional stiffness distribution 
(fig. 11(d)) higher than the minimum requirement. 

Figure 12 shows stiffness deviation versus cycle 
number for the three matching locations: the root 
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(fig. 12(a)), the point of taper initiation (fig. 12(b)), 
and the tip (fig. 12(c)). Early in the optimization 
process, the flapwise and torsional stiffnesses are both 
unmatched. After cycle 10, the matchings improve, 
and after 25 cycles, all three matchings are good. 
At the tip (fig. 12(c)), matching proves to be quite 
difficult. The torsional stiffness is the last to match. 
The reason for this is that the blade initial design is 
tapered, and it is difficult to place a thin-wall section 
in the space near and at the tip and still match the 
stiffness required on the upper level. 

Example 3: tapered planform (“ initially 

unmatched stiffnesses”). In the previous exam- 
ples, the starting points used matched stiffnesses. 
The purpose of the present example is to demon- 
strate how the IADS procedure behaves when it is 
started from an inconsistent set of stiffnesses (i.e., 
unmatched stiffnesses). The starting point for the 
optimization is shown in figure 8(c). The initial stiff- 
nesses used in the upper level are much larger than 
the stiffnesses obtained from the lower level design 
variables. The initial and final values for the blade 
planform, performance measures, and dynamics mea- 
sures are given in table 5(a). The initial and final 
constrained frequencies are included in table 5(b). 
The final upper and lower level stiffnesses are shown 
in figure 13, which also shows that the optimization 
procedure is able to match the upper and lower level 
stiffnesses successfully. Figure 14 shows the stiffness 
deviations for the three matching locations: the root 
(fig. 14(a)), the point of taper initiation (fig. 14(b)), 
and the tip (fig. 14(c)). As shown in the figure, 
the optimization procedure is able to match all three 
stiffnesses after 25 cycles, but it is at the expense of 
upper level performance. (See table 5.) The power 
required for all three flight conditions has increased 
substantially along with the hub shear. Notice that 
a bending frequency /j g has shifted frequency inter- 
vals. From these results it appears that although the 
optimization procedure will converge when starting 
from an initial point that has unmatched stiffnesses, 
starting with a set of consistent stiffnesses is better. 

Observations on Effect of Initial Design 
Study 

The IADS procedure has been exercised for three 
starting blade planforms: a rectangular planform 

with matched stiffnesses, a tapered planform with 
matched stiffnesses, and a tapered blade with un- 
matched stiffnesses. In all cases, the procedure is able 
to find converged feasible designs. When comparing 
examples 1 and 2 (tables 3 and 4, respectively), the 
reader will find two different final blade designs (i.e., 
the design variable values are different) with essen- 


tially the same objective function value. Apparently, 
many different combinations of design variables exist 
that satisfy the matching constraints, and more than 
one of these is optimal. The final solution depends 
on initial conditions. In example 3 (table 5), the 
optimizer appears to converge to a suboptimal solu- 
tion when compared with example 2. Both examples 
started from the same planform, but example 2 starts 
with matched stiffnesses and example 3 starts with 
unmatched stiffnesses. Because the initial matching 
of the stiffnesses is relatively easy, this suggests that 
the initial matching should always be enforced. 

When comparing all three examples, the reader 
will also notice that each initial blade has a different 
frequency range for the bending and torsional fre- 
quencies and that each final blade design has a fre- 
quency which has shifted a frequency interval (e.g., 
/ 6 6 in example 3). During the approximate analy- 
sis, the optimizer can change the upper level design 
variables such that a frequency can shift intervals. 
However, as the design- variable move limits are re- 
duced, this shifting is less likely to occur. 

At the present time, no automatic move-limit ad- 
justment is used in the upper level approximate anal- 
ysis. However, the stiffness deviation information 
(e.g., fig. 10) can be collected and used to determine 
when a move-limit adjustment is necessary for an up- 
per level design variable during the overall optimiza- 
tion process. 

Study on Effect of Coordination 

Parameter e 

The purpose of this study is to demonstrate the 
effect of e in the coordination constraint (eq. (26)) 
on the optimization procedure. Results for three 
values of e (0.4, -0.2, and -0.4) are presented in 
table 6 and figures 15, 16, and 9, respectively. If 
£ is a large positive value, the levels are essentially 
independent. The upper level is free to change 
the upper level stiffness and chord distributions in 
any way that will reduce the upper level objective 
function. The only requirement is that the overall 
stiffness matching should not degrade by more than 
the value of e from the best match found on the last 
lower level optimization. For example, if e = 0.4, 
the stiffness matching can degrade by 40 percent and 
still satisfy the coordination constraints. Therefore, 
the procedure could possibly converge with the upper 
and lower level stiffnesses being mismatched by as 
much as 40 percent. A negative value for £ means 
that the upper level must improve the matching 
achieved on the lower level by that amount. This 
section of the paper presents results for several values 
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of e using the starting point in figure 8(a) that is also 
used in example 1. 

One choice for e would be zero. This would mean 
that the upper level cannot degrade the matching 
achieved on the lower level. This value was found to 
be too restrictive for the optimization process, and 
the procedure converged in 3 cycles with very lit- 
tle change in the upper level design variables. The 
reason for this can be seen by examining the coordi- 
nation constraint (eq. (26)). At the start of the up- 
per level optimization, the coordination constraint at 
each matching location is active (i.e., g — 0) because 
F u is equal to As the upper level optimizer tries 
to change the upper level design variables, the coor- 
dination constraints become violated. Therefore, the 
upper level optimizer makes only small changes and 
the process converges in 3 cycles. 

As shown in table 6, when e = 0.4, the optimiza- 
tion process is able to improve the performance and 
dynamics measures over the initial blade values and 
improve the lower level (i.e., satisfy the stress con- 
straints). This improvement is achieved at the ex- 
pense of stiffness matching. Figure 15 shows the final 
stiffness distributions for the upper and lower levels. 
The lower level is not able to find a set of stiffnesses 
to match those required by the upper level. This fi- 
nal result is technically a feasible design because all 
the constraints are satisfied. Recall that the upper 
and lower stiffnesses need only be as close as possible 
(lower level objective function). The upper level co- 
ordination constraints do not require the upper and 
lower level stiffnesses to match exactly. 

When e = —0.2, the optimization procedure is 
able to obtain a design that has some improvement 
over the initial starting point (table 6). The upper 
level objective function is reduced slightly, but not 
as much as when £ is positive. As shown in figure 16, 
the upper and lower stiffnesses match well for the 
chordwise stiffness (fig. 16(a)), the flapwise stiffness 
(fig. 16(b)), and the torsional stiffness (fig. 16(c)). 
The lower level is able to obtain an extensional 
stiffness which is slightly larger than that required 
by the upper level. 

Of the values used in this work, the best value 
for e is —0.4 because the optimization procedure is 
able to obtain improvement on the upper level and 
find a set of consistent stiffnesses on the lower level. 
These results (example 1) are included in table 6 for 
completeness. The stiffness distributions are shown 
in figure 9. 


Observations on Effect of Coordination 

Parameter e Study 

As shown previously, positive values of £ result in 
upper level improvement, but poor stiffness matching 
and negative values of e result in both upper level 
improvement (although not quite as good as when £ is 
positive) and good stiffness matching. This suggests 
that a gradual reduction in £ from a positive value 
to a negative value could be beneficial. The IADS 
procedure was run with e = 0.4 for 8 cycles, £ = 0.2 
for 8 cycles, £ = -0.2 for 8 cycles, and e = -0.4 for 
8 cycles. This technique of gradually reducing the 
value of £ did not work. Presumably, the upper level 
planform area and upper level stiffnesses increased 
to improve the upper level objective function when £ 
was positive so that by the time that e was negative, 
the stiffness matching was achieved at the expense 
of performance and dynamic improvement on the 
upper level. This situation is analogous to example 3 
in which the mismatched initial conditions resulted 
in stiffness matching at the expense of upper level 
improvement. 

Study on Comparison of Two-Level and 

Single-Level Optimization Procedures 

The IADS procedure is compared with a more 
traditional optimization procedure without multi- 
level decomposition (i.e., the single-level optimiza- 
tion procedure). The single-level optimization pro- 
cedure combines local and global design variables 
and simultaneously evaluates aerodynamics, dynam- 
ics, and structures. The design variables are the same 
as those used in the IADS procedure with the excep- 
tion of the stiffness design variables ( EI XX , EI ZZ , GJ , 
and EA) which are no longer needed. These stiff- 
nesses are calculated from the wall thicknesses, the 
lumped areas, and the blade planform. The con- 
straints are the same constraints used in the IADS 
procedure with two exceptions. First, the coordi- 
nation constraints are no longer needed because the 
stiffnesses are calculated from the design variables. 
Second, the cumulative constraint, (KS), (eq. (25)) is 
no longer needed and the individual constraint com- 
ponents (eqs. (14)-(24)) are used. The single-level 
optimization procedure has 28 design variables and 
218 constraints (eqs. (2)— (11) and (14)-(24)). The 
optimizer consists of CONMIN and an approximate 
analysis. All derivatives are calculated by forward 
finite differences and results are presented for four 
cases. Cases 1 and 2 compare the single-level and 
multilevel approaches using the initial designs in fig- 
ures 8(a) and 8(b), respectively. Cases 3 and 4 in- 
vestigate whether the multilevel approach can im- 
prove on the best solutions found by the single-level 
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approach. All the IADS cases use a value of -0.4 for 
the coordination parameter (s). 

Case 1: rectangular starting design. The 
starting point for case 1 is shown in figure 8(a), and 
initial and final results for the two approaches are 
given in table 7. Both approaches show improvement 
over the initial design. The single-level approach 
finds a better overall design in terms of the objective 
function than the multilevel approach and is faster 
to converge. The multilevel approach has slightly 
less performance improvement but more dynamics 
improvement than the single-level approach. 

Case 2: tapered starting design. The start- 
ing point for case 2 is the tapered blade shown in 
figure 8(b) . Initial and final results for the two 
approaches are given in table 8. For this start- 
ing point, only the multilevel approach converges 
to a feasible design. The single-level approach was 
manually discontinued after 64 cycles. The final 
single-level design includes large lumped masses re- 
sulting in a weight constraint that is grossly violated 
(W = 4.4 lb). 

Case 3: starting from feasible single-level 
design. The purpose of case 3 is to see if the IADS 
optimization procedure can improve on the best so- 
lution obtained by a single-level optimization proce- 
dure. The initial design for the multilevel procedure 
is rectangular to 0.7131 percent radius and then ta- 
pers linearly to the tip with a taper ratio of 2.21 to 1. 
This is the best solution found using the single-level 
approach in case 1 (table 1). The results are given 
in table 9. The multilevel approach is able to im- 
prove the design only slightly over that obtained by 
the single-level approach. 

Case 4 : starting from infeasible single- 

level design. Recall in case 2 that the single-level 
approach did not find a feasible design although the 
IADS procedure did find a feasible design. The pur- 
pose of case 4 is to see if the IADS optimization pro- 
cedure can start from that final infeasible single-level 
solution (table 8) and obtain a feasible design. Thus, 
the initial design for the multilevel procedure has a 
pretwist of —13.22° and a planform that is rectan- 
gular to 0.4934 percent radius and then tapers lin- 
early to the tip with a taper ratio of 3.098 to 1 with 
a root chord of 0.4899 ft. The IADS procedure is 
able to find a feasible design, and the initial and fi- 
nal results are given in table 10. The final upper and 
lower level st iffne ss distributions for the three match- 
ing locations are shown in figure 17 and the stiffness 
deviations are shown in figure 18. 


Observations on multilevel versus single- 
level optimization cases. The IADS procedure is 
compared with a more traditional optimization pro- 
cedure without decomposition (single-level optimiza- 
tion procedure) for four cases. In the first and second 
cases, the initial design had a rectangular planform 
and a tapered planform, respectively. In the third 
and fourth cases, the IADS procedure is started with 
the best designs from a single-level optimization pro- 
cedure. The multilevel optimization approach is able 
to find feasible final designs regardless of the initial 
planform design. The single-level optimization pro- 
cedure could find only a feasible final design when the 
initial planform was rectangular (case 1). When the 
initial planform was tapered (case 2), the single-level 
optimization procedure could not find a feasible de- 
sign and was terminated manually because the blade 
weight constraint was grossly violated. In cases 3 
and 4, the multilevel approach was started from two 
single-level optimization results of cases 1 and 2, re- 
spectively. In case 3, the multilevel approach im- 
proved the single-level design slightly. In case 4, the 
initial design for the IADS procedure is the infeasible 
single-level design that has a grossly violated weight 
constraint. The IADS procedure is able to start with 
this highly infeasible design and find a feasible design. 

Concluding Remarks 

An integrated aerodynamic/dynamic/structural 
(IADS) optimization procedure for helicopter rotor 
blades has been developed. The procedure combines 
performance, dynamics, and structural analyses with 
a general-purpose optimizer using multilevel decom- 
position techniques. At the upper level, the structure 
interacts with the disciplines of aerodynamics and 
dynamics in terms of global quantities (stiffnesses 
and average strains). At the lower level, the struc- 
ture is defined in terms of local quantities (detailed 
dimensions of the blade structure and stresses). 

The IADS procedure consists of an upper level 
optimization, a lower level optimization, and a coor- 
dination task. The upper level objective function is a 
linear combination of performance and dynamic mea- 
sures. Upper level design variables include pretwist, 
point of taper initiation, taper ratio, root chord, 
blade stiffnesses, tuning masses, and tuning mass lo- 
cations. Upper level constraints consist of limits on 
power required in hover, forward flight, and maneu- 
ver; airfoil drag; minimum tip chord; trim; blade nat- 
ural frequencies; autorot at ional inertia; blade weight; 
and average strains. 

The lower level optimization sizes the internal 
blade structure to provide the stiffnesses required by 
the upper level and assure the structural integrity of 
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the blade. The lower level design variables are the 
box beam wall thicknesses and several lumped areas 
that are analogous to longitudinal stringers in a wing 
box cross section. The lower level objective function 
is a measure of the difference between the upper 
level stiffnesses and the stiffnesses computed from the 
wall thicknesses and lumped areas. The lower level 
constraints are on Von Mises stresses, extensional 
stiffnesses, thin-wall theory, and dimensional limits. 

The coordination task consists of a set of up- 
per level constraints that link the levels and promote 
consistency between the upper and lower level stiff- 
nesses. A coordination parameter is included in each 
constraint. This parameter specifies the amount of 
coupling between the levels. A proper value for the 
coordination parameter is found to be crucial to the 
success of the IADS procedure. If the parameter has 
a positive value, the procedure will converge but the 
final stiffness matching can be unacceptable. If the 
parameter value is too small (approximately zero), 
the optimization process will terminate without im- 
proving the dynamics or performance measures. A 
small negative value for the coordination parameter 
encourages the upper level to improve dynamics and 
performance using stiffness values that the lower level 
can match. 

The IADS procedure is demonstrated by using a 
model-size rotor blade for several initial blade plan- 
forms and varying amounts of coupling between the 
levels. In addition, the IADS multilevel procedure 
is compared with a more traditional optimization 


procedure without decomposition (a single-level op- 
timization procedure). In all cases, the IADS pro- 
cedure achieves successful results. It converges to 
a feasible design regardless of whether the initial 
design had a set of consistent stiffnesses. However, 
initializing the upper level stiffnesses with the stiff- 
nesses calculated from the lower level design variables 
greatly improves the final design. For the cases stud- 
ied, the IADS procedure is found to be superior to 
the single-level optimization procedure. The IADS 
procedure converges to a feasible design, even when 
the single-level procedure does not. Furthermore, the 
IADS procedure improves upon the best design found 
by the single-level optimization procedure. 

The IADS procedure exploits the couplings and 
interactions between the disciplines of aerodynamics, 
dynamics, and structures. It provides an efficient 
method to integrate structures and/or structural 
properties into an optimization procedure because it 
guarantees that a structure with a consistent set of 
structural properties can be found. The IADS proce- 
dure provides an optimization technique that is com- 
patible with industrial design practice in which the 
aerodynamic and dynamic design is performed at a 
global level and the structural design is carried out at 
a detailed level with considerable dialogue and com- 
promise among the groups. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
August 25, 1994 
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Appendix A 


(ref. 28) given by 


Coordination Constraint 

In a multilevel decomposition approach, the cou- 
pling between levels is done through a coordination 
procedure (e.g., refs. 20 and 25). In the present work, 
the coordination procedure based on reference 20 is 
used to reconcile the stiffnesses required on the upper 
level with the stiffnesses that the lower level can actu- 
ally obtain. This reconciliation results in one upper 
level constraint at each matching location 


dF L 

dF L 

T dKS 

dDY t 

n ” dDVi 

o 1 

n A dDVi 
0 * 


where dF ^ /dDVi is the derivative of the optimum 
lower level objective function with respect to the up- 
per level design variables, <9KS/<9DVi is the deriva- 
tive of the active lower level constraint (eq. (25)) with 
respect to the upper level design variables, X T is the 
Lagrange multiplier given by 


g = F u — (1 + s)F„ < 0 (Al) 

where Fq is the most recent value of the lower 
level objective function (i.e., the optimum value of 
eq. (13)), F u is an estimate of the change in Fq that 
would be caused by a change in the upper level design 
variable values, and e is the coordination parameter. 
This coordination parameter specifies how much the 
upper level can either degrade or improve the overall 
stiffness matching achieved on the lower level, and it 
may also be interpreted as a measure of how closely 
coupled the two levels are. If e has a positive value, 
the two levels are not closely coupled (i.e., they are 
essentially independent). The upper level can change 
the upper level stiffness and chord distributions in 
any way that will improve the upper level objective 
function as long as the stiffness matching is not 
degraded by more than the amount of e . If £ has 
a negative value, the two levels are closely coupled 
and the upper level is commanded to improve the 
matching by the amount of e. 

Equation (Al) is the general form of the coordi- 
nation constraint as formulated in reference 15. The 
form of the coordination constraint used in this work 
is obtained by approximating F u in terms of the cur- 
rent optimum lower level objective function (F 0 ). If 
Fq is expanded in terms of a first-order Taylor se- 
ries about the lower level optimum, then F ^ can be 
approximated by 


NDV 

f u = f 0 l +Y, 

i = 1 


dF L 

dDV t 


AD Vi 

0 


(A2) 


where 

dF L 

dDVi 
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DVj is an upper level design variable and 
is the total optimum sensitivity derivative 


where dKS/dv is the vector of derivatives of the ac- 
tive lower level constraints with respect to the lower 
level design variables at the lower level optimum. At 
a lower level optimum, will be positive, and if no 
lower level constraint is active, is set to zero. By 
substituting equation (A2) into equation (Al), the 
coordination constraint g is approximated by 


/ L ^ dF L 

9 ( F ° + S ^V. 

0 ADVi ) 

- (1 + e)F 0 L < o 

\ l ~ 1 

(A5) 

or simplifying gives 



/NDV , F L \ 

eFg < 0 (A6) 


From substituting equation (A3) into equation (A6), 
the coordination constraint becomes 


dF l 


dv 


(A4) 


NDV / q F L 

.t dKS 

\1 

[S ( 8dv < 

ADVl A dDVi 
0 1 

„)] 


-eFf < 0 


which is the form implemented in this work. 


(A7) 


The derivative of the coordination constraint is 
obtained by differentiating equation (A7) with re- 
spect to upper level design variables. Thus, 


dg _ dF L 

dDVi “ dDV t 


(A8) 
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Appendix B 
Drag Constraints 

Rotor blades operate over a wide range of flight 
conditions: hover, low-speed forward flight, high- 
speed forward flight, and maneuver. In addition, the 
blade encounters different conflicting phenoma as it 
rotates. Figure 19 shows a top view of the rotor disk 
with the advancing side when the blade is between 
4* = 0° and 180° and with the retreating side when 
the blade is between 4* = 180° and 360°. On the 
advancing side, the blade encounters a higher net 
velocity than it does on the retreating side where 
an area exists in which the flow is reversed. This 
reversed-flow region does not contribute any lift on 
the blade. In this work and in previous work (refs. 3, 
4, and 14), aerodynamic concerns such as drag di- 
vergence and blade stall are expressed in terms of 
constraints on the drag coefficient at various azimuth 
angles (4'). 

On the advancing side of the blade, c d is checked 
at every spanwise station. The constraint is formu- 
lated so that the largest c d is less than a given value of 
c d,al]- At a given azimuth angle, the largest section 
drag coefficient is selected from the aerody- 

namic stations along the blade span. (See fig. 20.) 
Thus, 


c d,m ax max ( c d 5 l » c d,2i • ■ • j c {f,MRA) 

($ = 0°, 15°, 30°,..., 180°) (Bl) 


where c di is the section drag coefficient in the zth 
aerodynamic segment and MRA is the total number 
of aerodynamic segments. 

Similarly, on the retreating side of the blade, c d 
is constrained at each 'P. The difference is caused by 
the reverse-flow region which occurs on the retreating 
side of the blade, and this must not be considered. 
In the reverse-flow region, c d is large because of the 
reversed flow. The velocities in this region have 
a tangential velocity similar to the advance ratio. 
Therefore, on the retreating side of the blade at a 
given azimuth angle, the largest drag coefficient is 
give by 


C flT,max max(q/£, c d,k+ 1> * • • » MRa) 

= 180°, 195°, 210°, . . . , 360°) (B2) 


where c d is the first drag coefficient corresponding 
to the first value of q ^ outside the reverse-flow region 
as shown in figure 21, 


At a given azimuth angle, the constraint is for- 
mulated as 


9 = 


d, max 
c d,all 


- 1 < 0 


(B3) 


where c d ^\ is the allowable drag coefficient and 
c d,max 15 S* ven by equations (Bl) or (B2). Because 
several airfoils may be used along the blade, a com- 
posite c d a ii is used. Different values of c d3 \i could 
be used for the advancing and retreating side con- 
straints, but in this work the same value of 0.12 is 
used. 


17 



Appendix C 

Lower Level Structural Analysis 

The purpose of this appendix is to summarize the 
elementary equations describing the geometric and 
structural analysis for the lower level structure. A 
typical cross section of the thin-walled isotropic box 
section is shown in figure 5. For simplicity, the top 
and bottom wall thicknesses, t\ and £3, respectively, 
are equal. The total cross-sectional area (A) is the 
sum of the cross-sectional areas of the box beam 
elements A % and the lumped areas aj (described in 
the main text). Thus, 

n m 

a = Aj + yy 0.7 (ci) 

i—l j = 1 

By using the familiar relations, the centroid of the 
cross section is calculated from the equations 


= 


n 


m 


y Ajx l + e a j%j 

i= 1 j = 1 


A 


(C2) 


and 


n m 


E Am 4- E *j*j 

i—l j = 1 



(C3) 


where x\ and z % are coordinates in the chordwise 
and flapwise directions, respectively, that specify the 
distance of the centroid of the zth element area {Ai) 
from the reference x - and 2-axes shown in figure 5. 
Similarly, Xj and Zj are coordinates that specify the 
distance of the centroid of the jth lumped area (aj) 
from the reference axes, n is the number of elements 
that the cross section is divided into for ease of 
calculations, and m is the number of lumped areas. 

Next, the area moments of inertia of each element 
about its centroidal x- and 2-axes are calculated from 


, h4 
x > k ~ 12 

(fc = 1,2, . . 

. , n + m) 

(C4) 

, h * b l 

*z,k - 12 

(k = 1,2,.. 

. , n + m) 

(C5) 


where b k is the base of the fcth rectangular element, 
h k is the height relative to the x-axis, and I XZj i — 0 
for symmetric elements. By using the parallel-axis 
theorem, the moments of inertia of each element are 


found with respect to the centroid of the box beam 
as 2 

Icx,k — -Tx,fc 1 (C6) 

hz,k = h,k + A k C 2 k J 

where and I cz ^ are the moments of inertia of 
the fcth element about the centroid of the box beam, 
I x k and I z k are the moments of inertia of the fcth 
element about its centroidal axes, and d k and c k are 
the distances from the centroid of the element to the 
centroid of the box beam in the x- and 2-directions, 
respectively. The total moments of inertia for the box 
beam are equal to the sum of the element inertias: 


I xx — ^ 

Izz = } ^ ^czM 


(C7) 


The polar moment of inertia (J) for the box beam 
is calculated by using the method described in refer- 
ence 38 which gives 


J = 


4 A 2 
§ ds/t 


(C8) 


where A c is the enclosed area of the mean periphery 
of the box beam wall, ds is the differential circumfer- 
encial length along the box beam, and t is the local 
thickness of the wall. 

In order to calculate the lower level objective 
function, the bending and torsional stiffnesses of the 
box beam are necessary. For an isotropic beam, 
the moments of inertia I X x an d I zz calculated in 
equation (C7) are multiplied by Young’s modulus E 
to acquire the bending stiffnesses EI X x an< 3 EI ZZ in 
the chordwise and flapwise directions, respectively. 
Similarly, the polar moment of inertia is multiplied 
by the torsional modulus of elasticity G to acquire 
the torsional stiffness of the beam GJ . 

The stresses for the constraints in the lower level 
optimization are evaluated at the corners of the box 
beam by using the Von Mises stress measure given 

by / 

V ( a , t) = Vct 2 + 3 t 2 (C9) 

where a is the axial bending stress at the outer fiber 
of the cross section, which is given by 


= (^f) Xouter + (1*7) “ outer + A (C10) 

and t is the shear stress due to torsion in the wall of 
the section with thickness t, which is given by 


2 ~Act 


(Cll) 
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where M zz is the flapwise moment, M xx is the lag 
moment, CF is the centrifugal force, and Mj is 
the torque at the section. The shear stress due to 
transverse loads has been neglected for simplicity. 


Here, M zz , M xx , N , and M? are computed in the 
upper level analysis for forward flight and maneuver, 
multiplied by a factor of safety (If), and then passed 
to the lower level. 
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Appendix D 
Wing Box Fitting 

At a given radial location, the outer dimensions 
of the box cross section are determined by placing a 
rectangular box of maximum area within the given 
airfoil cross section using a modified version of the 
method described in reference 34. As shown in 
figure 22, the outer dimensions (where b denotes 
the nondimensional box beam width and h denotes 
the nondimensional box beam height) depend on the 
airfoil section and the local chord c. (Note that the 
upper coordinates z u , lower coordinates z h horizontal 
coordinates x, and maximum thickness location t m ax, 
which are all normalized with respect to the chord c, 
are given in ref. 35 for the RC(4)-10 airfoil and in 
ref. 36 for the RC(3)-08 and RC(3)-10 airfoils.) 

The procedure for the wing box fitting is de- 
scribed below, and the nondimensional box beam 
height (h) is determined first. By starting from the 
leading edge of the airfoil at an initial point x/, winch 
is given as 


z £ and zj are determined by linearly interpolating 
between the respective upper and lower airfoil coor- 


dinates. The nondimensional box beam height ( h ) is 
given by 

h = z l u - z\ (D2) 

Next, by starting at a distance x r , where 

*r - 1 - ff») 

the location of the right side of the box is similarly 
determined. If the box is not within the airfoil shape, 
x r is increased from the trailing edge until the box 
is within the airfoil shape. The nondimensional box 
beam length (6) is given by 

b — x r — xi (D4) 

and the nondimensional area A of the box is given 

by 

A = bh (D5) 

Next, xi is incremented by 1 percent and the process 
is repeated to compute a new area. The larger of the 
two areas is kept. The process is repeated until three 
consecutive areas are within a given tolerance. When 
this occurs, a rectangular box of maximum area has 
been determined. 
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Table 1. Bounds for Design Variables Used in Optimization Examples 


Quantity 

Design variables 

Lower bound 

Upper bound 

Twist, deg 

-20.0 

-5.0 

Taper initiation, r/R 

0.26 

0.985 

Taper ratio 

0.05 

5.0 

Root chord, ft 

0.05 

0.833 

EI XX , lb-ft 2 

50.00 

20 000000.0 

EI ZZ , lb-ft 2 

5.00 

1000.0 

GJ, lb-ft 2 

5.00 

1000.0 

EA, lb 

1000.00 

200000000.0 

mi, slug/ft 

0 

0.50 

Vi,r/R 

0.24 

0.95 

t(, ft 

0.00008 

0.01 

ft 2 

0 

0.00004 




Table 2. Parameters and Flight Conditions Used in Optimization Examples 

(a) Parameters 


Minimum autorot at ional inertia, AI m i n , lb- ft 2 

Allowable drag coefficient, c^ a n 

Minimum tip chord, Ct >m: j n , ft 

Number of blades, N 

Number of aerodynamic segments: 

HOVT . . 

CAMRAD/JA 

Number of structural segments 

Number of design variables: 

Upper level 

Lower level 

Power available, P a , hp ......... 

Blade radius, P, ft 

Maximum blade mass, W, lb 

Factor of safety, If • 

Frequency increment, A/, per rev 

Allowable average strain, s a , ft/ft . . . . . 
ITER max * • * 

Pi nax ■ * •" • 7 - * 1 * 

Allowable stress, < 7 a , lb/ ft 2 

Young’s modulus, E, lb/ft 2 


23.69 
0.12 
0.083 
. 4 


19 

18 

50 


22 

18 (6 per location) 

20 

4.68 

3.5 

2.0 

0.1 

0.05 

40 

300 

, . . 8.352 x 10 6 

. . . 15.26 x 10 8 


(b) Flight conditions 

Rotational velocity (in Freon with density of 0.006 slug/ft 8 ), 17, rpm 
Hover tip Mach number 

C L : 

Hover 

Forward flight * 

Maneuver 

C D : 

Forward flight 

Maneuver 

Advance ratio: 

Forward flight 

Maneuver 


639.5 

0.628 


0.00810 

0.00810 

0.00985 


-0.000811 

-0.000596 


0.35 

0.30 
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Table 3. Example 1: Rectangular-Planform Starting Point With Matched Stiffnesses 


(a) Initial and final values for blade planform, performance measures, 
and dynamic measures 


Quantity 

Initial value 

Final value 

Hover power, hp 

14.81 

14.50 1 

Forward flight power, hp 

13.26 

12.96 

Maneuver power, hp 

12.22 

11.94 

Hub shear, lb 

2.1 

1.1 

Objective function . 

20.578 

19.9107 

Twist, deg 

-9.0 

-11.47 

Taper initiation, r/R 

0.7 

0.7010 

Taper ratio 

1.0 

1.664 

Root chord, ft 

0.3449 

0.3770 

mi, slug/ft 

0 

0.00027607 

m 2 , slug/ft 

0 

0.0031988 

m 3 , slug/ft 

0 

0.0020144 

VI, r/R 


0.4503 

S /2 ,r/R 


0.5830 

S/3 ,r/R 


0.4534 

Cycles to converge 


76 


(b) Initial and final values for constrained frequencies 


Frequency, per rev 

Initial value 

Final value 

fb,3 

2.60 

2.68 

fb, 4 

3.77 

4.57 

fb, 5 

4.52 

4.88 

fb, 6 

7.22 

7.55 

ft, 1 

7.30 

7.30 

ft, 2 

3.61 

3.83 


(c) Final values for lower level design variables and upper level stiffnesses 


Variable 

Radial location 1 
(root) 

Radial location 2 
(point of taper initiation) 

Radial location 3 
(tip) 

Final lower level design variables 

<i,ft 

0.002366 

0.002427 

0.0004517 

t 2 , ft 

0.003261 

0.009954 

0.0003766 

£ 4 , ft 

0.003414 

0.009954 

0.0003766 

ai, ft 2 

0.00003341 

0.00003293 

0.00001610 

^ 2 , ft 2 

0.00001615 

0.00003084 

0.00001192 

03 , ft 2 

0.00003281 

0.00003293 

0.00001610 


Final upper level stiffnesses 


El xx , lb-ft 2 . . . 

2057.0 

2974.1 

153.73 

EI ZZ , lb-ft 2 . . . 

122.21 

140.03 

8.6135 

GJ, lb-ft 2 .... 

127.93 

128.53 

5.8743 

EA, lb 

797370 

1647300 

212 230 
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Table 4. Example 2: Tapered-Planform Starting Point With Matched Stiffnesses 


(a) Initial and final values for blade planform, performance measures, 

and dynamic measures 


Quantity 

Initial value 

Final value 

Hover power, hp 

14.85 

14.74 

Forward flight power, hp 

13.38 

13.02 

Maneuver power, hp ........... - 

11.93 

11.84 

Hub shear, lb 

0.6 

0.66 

Objective function 

19.876 

19.9326 

Twist, deg . 

-16.0 

-10.85 

Taper initiation, r/R 

0.8 

0.37 

Taper ratio 

3.0 

1.636 

Root chord, ft 

0.45 

0.4932 

mi, slug/ft 

0 

0.008961 

7712 , Slug/ft 

0 

0.01354 

m3, slug/ft 

0 

0.0246 

yi,r/R ■ • 


0.24 

y 2 ,r/R . . . 


0.6164 

2/3 ,r/R 


0.6215 

Cycles to converge 


93 


(b) Initial and final values for constrained frequencies 


Frequency, per rev 

fb, 3 

fb, 4 

fb, 5 

fbfi 

ft, l 

ft, 2 


Initial value 


2.93 

5.64 

6.22 

10.25 

7.30 

6.45 


Final value 


2.86 

5.33 

6.68 

9.16 

7.30 

6.12 




Table 5. Example 3: Tapered-Planform Starting Point With Unmatched Stiffnesses 

(a) Initial and final values for blade planform, performance measures, 
and dynamic measures 


Quantity 

Initial value 

Final value 

Hover power, hp 

14.85 

16.64 

Forward flight power, hp 

13.27 

17.46 

Maneuver power, hp 

11.89 

14.89 

Hub shear, lb 

0.186 

2.45 

Objective function 

20.005 

24.624 

Twist, deg 

-16.0 

-11.98 

Taper initiation, r/R 

0.8 

0.8893 

Taper ratio 

3.0 

1.3148 

Root chord, ft 

0.45 

0.7364 

mi, slug/ft 


0.008546 

7712, Slug/ft 


0.0077966 

m 3 , slug/ft 


0.0090299 

yi,r/R 


0.32257 

V 2 ,r/R 


0.43886 

2/3 ,r/R 


0.39256 

Cycles to converge 


92 


(b) Initial and final values for constrained frequencies 


Frequency, per rev 

Initial value 

Final value 

fb, 3 

2.87 

2.90 

fb, 4 

5.54 

5.87 

fb, 5 

8.62 

8.10 

fb,6 

9.65 

10.50 

/t,l 

7.30 

7.30 

ft, 2 

5.48 

5.12 





Table 6. Effect of Coordination Parameter (e) on Multilevel Optimization Procedure 


Quantity 

Initial value 

Final value at e of — 

0.4 

-0.2 

-0.4 

Hover power, hp 

14.81 

14.44 

14.60 

14.50 

Forward flight power, hp . . . 

13.26 

12.77 

13.11 

12.96 

Maneuver power, hp 

12.22 

11.75 

11.96 

11.94 

Hub shear, lb 

2.1 

0.2072 

1.85 

1.1 

Objective function 

20.58 

19.48 

20.22 

19.9107 

Twist, deg 

-9.0 

-13.32 

-11.12 

-11.47 

Taper initiation, r/R .... 

0.7 

0.7859 

0.8246 

0.7010 

Taper ratio 

1.0 

3.155 

1.410 

1.664 

Root chord, ft 

0.3449 

0.3651 

0.3606 

0.3770 

mi, slug/ft 

0 

0.02571 

0.00135 

0.0002761 

m 2 , slug/ft 

0 

0.00211 

0.0000995 

0.0031988 

m 3 , slug/ ft 

0 

0.00099 

0.0000727 

0.0020144 

Vl,r/R 


0.4124 

0.3115 

0.4503 

yo.r/R 


0.4154 

0.3950 

0.5830 

yz,r/R 


0.4382 

0.4292 

0.4533 

Cycles to converge 


90 

152 

76 
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Table 7. Case 1 : Comparison of Single-Level and Two-Level Optimization 
Procedures for Rectangular Starting Design 


Quantity 

Initial value 

Final value for approach— 

Single level 

Two levels 

Hover power, hp 

14.81 

14.42 

14.50 

Forward flight power, hp . . . 

13.26 

12.87 

12.96 

Maneuver power, hp 

12.22 

11.83 

11.94 

Hub shear, lb 

2.1 

1.27 

1.1 

Objective function 

20.578 

19.5469 

19.9107 

Twist, deg 

-9.0 

-13.25 

-11.47 

Taper initiation, r/R .... 

0.7 

0.7131 

0.7010 

Taper ratio 

1.0 

2.21 

1.664 

Root chord, ft 

0.3449 

0.3813 

0.3770 

mi, slug/ft 

0 

0.003065 

0.00027607 

m 2 , slug/ft 

0 

0.04554 

0.0031988 

m 3 , slug/ft 

0 

0.013545 

0.0020144 

yi,r/R 


0.3073 

0.4503 

2/2 ,r/R 


0.5132 

0.5830 

2 / 3 , r/R 


0.5061 

0.4534 

Feasible design 


Yes 

Yes 1 


Table 8 . Case 2 : Comparison of Single-Level and Two-Level Optimization 
Procedures for Tapered Starting Design 


Quantity 

Initial value 

Final value for approach — 

Single level 

Two levels 

Hover power, hp 

14.85 

14.41 

14.74 

Forward flight power, hp . . . 

13.38 

12.67 

13.02 

Maneuver power, hp 

11.93 

11.68 

11.84 

Hub shear, lb 

0.6 

0.7137 

0.66 

Objective function 

19.876 

19.377 

19.9326 

Twist, deg 

-16.0 

-13.22 

-10.85 

Taper initiation, r/R .... 

0.8 

0.4934 

0.37 

Taper ratio 

3.0 

3.098 

1.636 

Root chord, ft 

0.45 

0.4899 

0.4932 

mi, slug/ft 

0 

0.002305 

0.008961 

m 2 , slug/ft 

0 

0.1636 

0.01354 

m 3 , slug/ft 

0 

0.2922 

0.0246 

Vl,r/R 


0.2420 

0.24 

VI ,r/R 


0.2420 

0.6164 

V3,r/R 


0.3610 

0.6215 

Feasible design 


No 

Yes 





Table 9. Case 3: Multilevel Approach Starting Prom Feasible 
Single-Level Optimization Solution 


Quantity 

Initial value 
from single- 
level approach 

Final value 
from two- 
level approach 

Hover power, hp 

14.42 

14.44 

Forward flight power, hp . . . 

12.87 

12.89 

Maneuver power, hp 

11.83 

11.83 

Hub shear, lb 

1.27 

1.13 

Objective function 

19.5469 

19.814 

Twist, deg . 

-13.25 

-12.51 

Taper initiation, r/R .... 

0.7131 

0.7450 

Taper ratio 

2.21 

2.4041 

Root chord, ft 

0.3813 

0.3834 

mi, slug/ft 

.003065 

0.001920 

m 2 , slug/ft 

0.04554 

0.08434 

m 3 , slug/ft 

0.013545 

0.01985 

y\,r/R 

0.3073 

0.3487 

1 12, r/R 

0.5132 

0.4397 

VZ,r/R 

0.5061 

0.4421 

Feasible design 

Yes 

Yes 


Table 10. Case 4: Multilevel Approach Starting Prom Infeasible 
Single-Level Optimization Solution 


Quantity 

Initial value 
from single- 
level approach 

Final value 
from two- 
level approach 

Hover power, hp 

14.41 

14.83 

Forward flight power, hp . . . 

12.67 

13.35 

Maneuver power, hp 

11.68 

11.96 

Hub shear, lb 

0.7137 

1.29 

Objective function 

19.377 

20.328 

Twist, deg 

-13.22 

-10.80 

Taper initiation, r/R .... 

0.4934 

0.36 

Taper ratio 

3.098 

2.164 

Root chord, ft 

0.4899 

0.5926 

mi, slug/ft 

0.002305 

0.0009530 

m 2 , slug/ft 

0.1636 

0.04803 

m 3 , slug/ft 

0.2922 

0.1042 

y\ ,r/R 

0.2420 

0.2913 

V 2 ,r/R 

0.2420 

0.3221 

V3 ,r/R 

0.3610 

0.4373 

Feasible design 

No 

Yes 




Upper level analysis 



Figure 1. General multilevel optimization procedure with two levels. 


















Figure 3. Flowchart of upper level analysis. 
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(a) Example 1: rectangular planform with matched upper and lower level stiffnesses. 




(c) Example 3: tapered planform with unmatched upper and lower level stiffnesses. 

Figure 8. Starting points. 
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(b) Point of taper initiation. 
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Figure 10. Convergence history of upper and lower level stiffness deviations for example 1. 
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Figure 1 1 . Final upper and lower level stiffness matching along blade radius for example 2. 






Figure 12. Convergence history of upper and lower level stiffness deviations for example 2. 
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Figure 14. Convergence history of upper and lower level stiffness deviations for example 3. 
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Figure 15. Final upper and lower level stiffness matching along blade radius for £ = 0.4. 
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Figure 16. Final upper and lower level stiffness matching along blade radius for e = -0.2. 
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Figure 18. Convergence history of upper and lower level stiffness deviations for case 4. 
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Retreating side 


Advancing side 



Figure 19. Rotor disk flow field. 
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